Exact Chiral Invariance at Finite Density on Lattice 
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Abstract. A new lattice action is proposed for the overlap Dirac matrix with nonzero chemical 
potential. It is shown to preserve the full chiral invariance for all values of lattice spacing exactly. 
O | It is further demonstrated to arise in the domain wall formalism by coupling the chemical 

D , potential count only to the physically relevant wall modes. 



1. Introduction 

A fundamental aspect of the phase diagram of Quantum Chromo Dynamics (QCD) in the T-\xb 
plane is the critical point, where T and \ib denote the the temperature and baryonic chemical 
potential respectively. Theoretical [IJ [2j [3] as well as experimental searches for locating it are 
currently going on [3]. Its discovery would be exciting in many ways. Apart from becoming a 
new milestone in our understanding of the nature of the strongly interacting matter, it would 
also be unique compared to the other known phase diagrams in the way theory and experiment 
compliment each other in locating the critical point in it. It is desirable to predict the existence 
of the QCD critical point (or rule it out) starting from first principles, employing only the QCD 
Hamiltonian. Lattice QCD is the best suited tool for such an exploratory exercise in view of its 
success in other non-perturbative aspects. 

The prime reason for expecting the critical point are the chiral symmetries QCD has due 
to the quark mass spectra. Our world has two light quarks, having masses much smaller than 
AqcDj an d a moderately heavy quark, making it a '2+1' flavour QCD. All attempts to explore 
nonzero /ib on the lattice use the staggered quarks. They have a U(l)y x U(1)a chiral symmetry 
on the lattice, and a corresponding order parameter, the chiral condensate. Staggered quarks, 
however, break the flavour and spin symmetry on the lattice, and have no flavour singlet axial 
anomaly. In the continuum limit of a — >■ 0, they have 4 quark 'flavours' of the same mass. 
The '2+1' flavour QCD is approximated by taking square root and fourth root of the staggered 
quark determinant. Whether this can restore the correct flavour symmetry as well as the chiral 
anomaly in the continuum limit has been a controversial issue. Moreover, it has been argued 
the 'rooting problem' becomes even worse [5] at nonzero hb- Since the QCD critical point is 
expected [6] to exist if only if one has two light flavours, and the flavour singlet anomaly is mildly 
temperature dependent [7], as in our world, it appears desirable to improve upon them. The 



overlap Dirac fermions |§J, or the closely related domain wall fermions [9], offer such a possibility 
to improve. Indeed, the overlap quarks have all the symmetries of the continuum QCD and also 
have an index theorem [TO] as well, raising the hope that even the anomaly effects could be 
well treated. Unfortunately, adding the chemical potential turns out to be nontrivial for them. 
Bloch and Wettig [11] made a proposal to do so but it violates [12] the exact chiral invariance 
on lattice as does the simple addition [13] of a baryon number term. In this talk I report on our 
alternative action |14] which does have exact and continuum-like chiral invariance on lattice for 
any value of the lattice spacing and any chemical potential. Its chiral order parameter permits, 
in principle, the task of mapping out the T-fis phase diagram, assuming that the algorithmic 
developments can handle the fermion sign problem well. 

2. Exact Chiral Invariance for \i ^ 

We begin by noting that the massless continuum QCD action can be written in a form where 
the chiral symmetry is manifest in terms of the fields appearing in the action and ask [TJ] if this 
can be imitated on lattice. The continuum QCD action at finite density in an explicit chiral 
invariant form is 

Sqcd = J d 3 x dr$ L (p + /j 7 4 )V>l + MP + M7 4 Wfl " , (1) 

where tp L = (1 - Ts)^/ 2 and ip R = (1 + Ts)^/ 2 witn $L = ^(1 + 7s)/ 2 and ip R = tp(l - 7s)/ 2 - 
The second term in the action is for gluons and will play no role in our discussion below. We 
assume that some usual convenient lattice form for it has been chosen and focus on the lattice 
fermionic action only. We propose that addition of chemical potential on the lattice be done 
in such an explicit chiral symmetry preserving manner as well, opting therefore for the overlap 
quarks on the lattice. 

The overlap quarks have all the symmetries of the continuum QCD but also have a nonlocal 
action: 

Sf = ^2$n aDov,nm tpm , ( 2 ) 
n,m 

where the sum over n and m runs over all the space-time lattice sites, a is the lattice spacing, 
and the overlap Dirac matrix D ov is defined by aD ov = 1 + 75 sgn (75 Dy/). sgn denotes the sign 
function. Dy/ is the standard Wilson-Dirac matrix on the lattice but with a negative mass term 
M G (0,2): 

aD w (x,y) = (4 - M)8 x , y - £[£/t(x - i)5 x _ ly ±±^ + l^Hu^S^J . (3) 

i=l 

The overlap Dirac matrix satisfies Ginsparg- Wilson relation [15], {75, D} = aDj^D and has 
exact chiral symmetry on lattice. The corresponding infinitesimal chiral transformations [16] 

Sip = £0:75(1 — aD ov )ip and dtp = ia"0 7 5 . (4) 

The generators of the transformation in Eq.(j2|) satisfy 75 = 1 and 75 = [75(1 — aD ov )] 2 = 1. 
One therefore defines the left-right projections for quark fields as ^ = (1 — 75 )^/ 2 and 
ipR = (1 + ■y$)ip/2, leaving the antiquark field decomposition as in the continuum. Such a 
decomposition is commonly done for writing chiral gauge theories [T7] on the lattice and is 
possible since tp and ip are independent fields in the Euclidean field theory. In analogy with the 
Eq.(Tf]) of continuum QCD, the action for the overlap quarks in presence of nonzero chemical 



potential may now be written down as 

S = ^2[i , n,L( aD ov + «M7 4 )nm'0m,-L + 4>n,R( aD ov + a^ A ) nm 1p m ^ R } (5) 
n,m 

= ^2i>n[(l ~ a^/2)aD ov + a^ 4 ] nm ip m . (6) 

n,m 

It is easy to verify that this action is invariant under the chiral transformation Eq. (J3|) for all values 
of afi and a and it reproduces the continuum action in the limit of a — > with fx — > \xjM scaling. 
In order to obtain the order parameter for checking if the symmetry is spontaneously broken, one 
usually adds a linear symmetry breaking term. Adding a quark mass term, am(ij)RipL + ipL^R,), 
one obtains the order parameter valid for all T and \i on the lattice by taking a derivative of 
the log of the partition function with respect to am as, 

Oijjifi) = lim lim (Tr — — j) 

am->o v->oo aK ov + am + a/i7 4 

where K ov = D ov (l — aD ov /2)~ 1 , such that {75, -fQ™} = 0. Although the discussion above is 
for a single flavour of quark, i.e, U(1)l x U(1)r symmetry, its generalization to Nf flavours is 
straightforward. Indeed, since it relies only on the spin-structure, the flavour index as well as 
the corresponding generator matrices just carry through. 



(7) 



3. Physical Picture 

Domain wall fermions are known to be akin to the the overlap fermions, and may provide a 
physical picture for introduction of fi in the way above. Moreover, full QCD simulations with 
them are easier, so it is also practically useful to check how the above exact chiral invariance 
can be preserved. The domain wall action of [IS] for fi = 0, is 

N S 

S = a V(z,s) [am 5 XfX > (d s ,i6 s ',N 5 P+ + P-Ss,n 5 $s',i) (8) 

x,x' s,s'=l 

-&x,x' (P-5 s ',s+l + P+Sa',8-1) + {a5P>w{x, x') + 8 XjX ') 5 S>S >] ip(x', s') , 

where P± = (1 ± 7s)/2 and N$, 05 are the number of sites and the lattice spacing in the fifth 
direction respectively. The physically relevant 4D fermion field is identified with the fermion 
fields at the boundaries of the fifth dimension as, 

= P-ip(x,l) +P + i(j(x,N 5 ) , t(ar) = ^(x, 1)P+ +^(x,iV 5 )P_ . (9) 

These are the lowest energy field configurations of the five dimensional theory and these 
correspond to the four dimensional chiral fermions in the limit m — > and ./V5 — > 00. In order 
to get a four dimensional fermion determinant representing the dynamics of the physical quark 
fields, one has to integrate out all the fermion degrees of freedom from the five dimensional 
action and remove the contribution of the higher energy configurations. The latter are the 
bulk modes since these are delocalized in the fifth dimension even when ./V5 — > 00. Integrating 
out the fermion degrees of freedom is conveniently done by recasting the domain wall action 
in terms of the fields m localized on four dimensional branes existing at each site i along the 
fifth dimension [TS]. These are related to the 5D fermion fields in an asymmetric way. The 
neighbouring pairs of rj are related by a transfer matrix, T = (1 + a^,HwP+)~ 1 (l — a^HwP-) 
with H\y = j^Dyy. Integrating over the 77 fields successively one obtains a determinant which 
is a functional of this transfer matrix: 

D^\ma) = det \p_ - maP+ - T~ Nb (P+ - maP_)l . (10) 



It has the contribution of all the fermions in the five dimensional theory. Choosing ma = 1, 
one obtains contribution from all but the lowest energy configurations which is divided out to 
obtain physical quark determinant. Taking first a§ — >■ limit and then N§ — > 00, the ratio of 
determinants turns out to be the determinant of the traditional overlap Dirac matrix. 

To get correct physics at finite /i one needs to construct the number density corresponding 
to the physical degrees of freedom. We proposed [2] that the following term be added to the 
domain wall action in Eq. (JS]) for nonzero density: 

N 5 

a 4 tp(x,s)afj, 74 S X)X i (Ss^Sgt^P- + P + 5 StN5 5 s/tN5 ) ip(x',s'). 

x,x' s,s'=l 

We then follow the same procedure outlined above for the fj, = case and rewrite the full action 
in terms of the 77-fields. The //-dependent term then becomes as below: 

\ia fji(a 5 H w P- - l) _1 74-P-7?i - fj N5 (a 5 H w P^ - l)' 1 -f i P + r] 1 (11) 

It is important to note that unlike the Bloch-Wettig case [20], the chemical potential in our 
proposal is coupled only to the physical fermions which are strongly localized on the boundaries 
of the fifth dimension, hence our transfer matrix continues to remain the same, i.e., is fi- 
independent. Moreover since we did not couple \i to the number density of the bulk modes, 
the pseudo-fermion action remains the same as in the \i = case. Integrating over the 77 fields 
one obtains the physical fermion determinant in a ratio form similar to that for fj, = with the 
det (ma) in the numerator generalised to, 

det D^ 5 \ma, /j,a) = det P_ — maP + + a\i (a^Hy/P- — l) -1 74 -P- 

- T~ Ns (P+ - maP. - an (a r> H w P + + l)" 1 74P+)] . (12) 

and the pseudo-fermion determinant in the denominator remaining unchanged. On taking the 
limits 05 — > and N5 — > 00, the ratio of determinant simplifies to det[D ov + (1 — D ov /2)(ma + 
a/ry 4 )], where both the dimensional parameters /i and m have been scaled by a factor of 1/M 
due to the tree level wavefunction renormalization [21]. A little algebra shows that 7 4 can be 
commuted through in the determinant above to yield the same overlap matrix of Eq. (|6|) with 
exact chiral symmetry on the lattice. 

The action in Eq.([6]) leads to an overlap fermion determinant which is identical to that in the 
recent work |22] with fermionic sources in the overlap formalism of |23j. The main difference 
is, however, in the necessity of sources in [22 to define the chiral symmetry. Indeed, the chiral 
symmetry transformation there is local, defined as rotation of the sources, while our Eq.Q 
is nonlocal, defined as the rotation of quark fields. The left-right symmetry is in-built in the 
formalism there whereas we needed to introduce the left-right projections in form of L- and 
P-fields to do so. 

It is easy to verify that our D ov (afj,), defined above, is not 75-hermitian. In general, therefore, 
it is not clear whether it may be diagonalizable. Noting that an M-scaling was essential in the 
continuum limit of our D ov with o/i/M governing the density dependent term, it is easily seen 
that for any arbitrary value of fi the density term can be made small by choosing small enough 
a. Hence it suffices to study the chiral anomaly at 0((afi/M) 2 ). The leading term in such 
an expansion in a/i/M, is D ov (0) which is 75-hermitian. Its eigenvalues come in complex pairs 
(A, A*), with the corresponding eigenvectors related by a 75-rotation. Using these properties, 
we showed that chiral anomaly arises as usual, and is governed by the number of zero modes of 
D ov(fl) even on the lattice. In particular, the contribution of the eigenvectors for the eigenvalue 



two, i.e, the 'doublers', to chiral anomaly is identically zero while that of the pairs (A, A*) cancels 
mutually. Using first order perturbation theory in a/i/M, it is easy to show that that the latter 
statement remain unchanged even for nonzero pna fine enough lattice for our proposal while 
the former is true on any afx/M. Since zero modes continue to be eigenvalues of 75, all one 
needs to do is to compute the trace of 75 in a new basis, rotated by a unitary transformation. 
After some algebra, one can show that the chiral anomaly remains unchanged on lattice for 
arbitrary fi for our proposal to introduce chemical potential for the overlap quarks, exactly as 
in the continuum. 

4. Summary 

In conclusion, we used the analogy with continuum QCD to demand a manifest chiral invariance 
in terms of L and .R-fields on the lattice to obtain aD ov + a//7 4 (l — aD^/l) as the exact chiral 
invariant form of the overlap Dirac matrix for any \i and/or lattice spacing. We wrote down the 
exact chiral order parameter on the lattice for chiral symmetry restoration in the entire T — jjlb 
plane. Chiral anomaly remains unaffected for our proposal for small enough lattice spacing 
a, just as in continuum. We provided a physically appealing interpretation of our proposal 
by explicitly writing down the corresponding domain wall fermion action at nonzero chemical 
potential and taking the limit of infinite fifth dimension. 
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